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1,  Introduction 

The  solution  of  the  problem  of  the  scattering  of  radio  waves  by  aggre- 
gates of  discrete  scatterers  in  random  motion  ('chaff,  rainfall,  etc.)  is  by  now 
classical^-'*'-].  The  time  structure  or  fading  characteristics  of  the  scattered 
radiation  is  caused  by  Coppler  shifting  of  the  incident  radiation  by  the  separate 
scatterers,  and  is  thus  simply  related  to  the  velocity  distribution  of  the  scat- 
terers. In  recent  years,  the  model  of  discrete  scatterers  has  been  borrovred  to 
explain  the  fading  of  radio  waves  scattered  by  dielectric  turbulence.  This  has 
led  to  the  models  of  Booker  and  Gordon L-',  Rice'-  ■' ,  and  others,  in  which  the  tur- 
bulence, at  least  for  the  piirposes  of  fading  calculations,  is  represented  as  an 
aggregate  of  'scattering  blobs',  with  velocities  resulting  from  turbulent  fluctua- 
tions and  convection  by  the  mean  wind.  Stich  models  lead  immediately  to  the  linear 
relation  between  fading  rate  and  carrier  frequency  which  is  characteristic  of 
Doppler  shifting.   In  this  paper  we  show  how  this  model  is  modified  when  the  speci- 
fically turbulent  inner  structure  of  the  random  dielectric  constant  field  is  taken 
into  account,  in  particular  when  the  dielectric  turbulence  obeys  the  Kolmogoroff- 
Obukhoff  structure  law  as  described  elsewhere'--'.  It  is  found  that  the  dependence 
of  fading  rate  on  carrier  frequency  differs  from  that  given  by  the  scattering  blob 
models. 


Suppose  an  antenna  A  emits  monochromatic  electromagnetic  radiation  of 
wavelength  X,  which  is  scattered  into  a  receiving  antenna  A  ,  at  a  distance  D 
beyond  the  geometrical  horizon  from  A,  by  the  dielectric  noise  in  the  common 
volume  V  of  the  antenna  beams.  For  simplicity,  and  without  loss  of  generality, 
we  assume  that  V  is  a  cube  centered  at  the  origin  0  of  a  cartesiaii  coordinate 
system  which  is  fixed  with  respect  to  the  earth.  Let  n  be  a  unit  vector  along 
the  line  AO  and  n'  a  unit  vector  along  the  line  OA  .  Then  the  scattering  wave 
vector  K  is  the  vector 

K  =  (2nA)  (n'-  n)  . 

In  this  section  we  assemble  preliminary  material  for  the  considerations  that 
follow;  a  more  detailed  discussion  of  this  material  with  the  same  notation  is 
given  by  the  author  ^  -' . 
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In  terms  of  the  scattering;  angle  ©  (see  Figure  l),  we  have* 

K  =   |k|   "  (hn/X)  sin(0/2)  ^   2ne/\     =  2nef/c  , 

where  f  is  the  carrier  frequency,  and  c  is  the  velocity  of  light.   (We  neglect 
the  variation  of  the  antenna  propagation  factors  in  V,  by  setting  them  equal  to 
xinity  everyi-jhere  in  V.) 


Figure  1 
Geometrical  configuration  of  a  radio  scattering  experiment  (schematic). 

The  received  scattered  radiation  will  be  described  by  giving  the  statis- 
tical structure  of  the  random  function  I(t),  the  complex  amplitude  of  the  current 
in  the  receiving  antenna  load  measvired  in  units  of  peak  amplitude  of  the  cior- 
rent  that  would  appear  in  the  receiving  antenna  load  if  the  antennas  were  in 
free  space  a  distance  D  apart  and  oriented  with  their  directions  of  maximum 
gain  along  the  line  joining  them.  We  recall  the  following  facts  about  l(t): 

We  shall  -use  the  symbol  '— '  to  mean  approximately  equal  to,  sometimes  only  to 
within  an  order  of  magnitude. 


1.  To  a  good  approximation,  l(t)  is  a  complex  gaussian  random  function 
of  time,  vrith  mean  zero.   Accordinr^ly,  I(t)  is  specified  (at  least  to  the  extei.t 
that  it  is  gaussian)  by  giving  its  covariance  <^l(t)I  (t+t")^  .   (The  asterisk 
denotes  the  complex  conjugate,  and  the  angular  parentheses  mean  the  ensemble 
average.) 

2.  I(t)  is  approximately  stationary  for  time  intervals  of  several 
minutes  duration.  Thus  the  covariance  of  l(t)  is  s\ibstantially  independent  of 
t  for  minutes  at  a  time,  i.e., 

<i(t)i*(t*r)>  --  <|i|^>   c(r)  , 

where  the  correlation  function  0(2^)  is  a  function  only  of  the  difference  variable 
r  =  t-t'. 

By  the  dielectric  noise  /\  e(r,t)  we  mean  the  random  departure  from  its 
mean  of  the  relative  dielectric  constant  at  the  point  (r,t).  We  make  the  follow- 
ing assumptions  about  ^  e(r,t): 

3.  /\^  e(r,t)  is  approximately  stationary  in  space  within  V  and  in  time 
•for  intervals  of  the  order  of  several  minutes  (as  is  l(t));  thiis  we  can  write  the 

covariance  of  ^  e(r,t)  as 

where  the  correlation  function  R(p,?^)  is  a  function  only  of  the  coordinate  dif- 
ferences  p  =  r  -  r  and  T  =  t  -  t  . 

U.  For  every  fixed  t,  the  random  fiinction  /s^  e(r,t)  is  substantially 
auto-correlated  only  over  a  distance  L  ,  the  macro-eddy  dimension  of  turbvilence 
theory  |  it  is  difficult  to  define  L  to  witldn  better  than  an  order  of  magnitude. 

Two  other  ass-umptions  that  we  shall  make  are: 

5.  The  magnitude  of  the  scattering  vector  satisfies  the  inequality 
KL  »  2n.  In  qualitative  language,  the  eddies  of  wave  number  K,  the  so-called 
scattering  eddies,  are  at  least  an  order  of  magnitude  smaller  than  the  macro- 
eddies.  This  condition  is  needed  if  statistical  turbxilence  theory  is  to  be  used 
to  make  deductions  about  ^^^  e(r,t)  and  consequently  about  l(t). 

6.  The  underlying  turbulent  velocity  field  v(r,t),  wlilch  convects  the 
dielectric  noise,  has  the  saj.ie  stationarity  properties  and  correlation  distance  as 


We  assume  some  acquaintance  on  the  part  of  the  reader  with  the  terminology  and 
lulence 

n-t...M 


results  of  statistical  turbulence  theory.  Basic  references  are  G.K,  Batche^or •- -■ 


and  L.  Landau  and  E.  Liftshitz' 


/\  6(r,t).  (Thus  must  be  the  case  if  the  dielectric  noise  is  caused  by  turbtdent 
mixing'--'-' ,) 

The  space-time  spectral  density  Y(k,co)  of  the  random  process  /\  e(?,t) 
and  the  space-time  correlation  function  R(p,T)  are  related  by  the  familiar  Wiener- 
Khintchine  formulas 

Y(k,co)   =  (2n)"^   j  R(p,r)  exp[i(k.p+wr)]  d?  dT  , 
R(p,'^')   =   )  Y(k,a))  exp[.i{k«p  +  ur)]  d^  doj  , 

where  all  variables  of  integration  range  from  -co  to  oo.  Note  that  Y(k,a3)  is 
normalized,  i.e.,  that 

I  Y(k,a))  dk  dw  =  1  , 
Of  particular  interest  is  the  qviantity  ^  (k,T)  defined  as 

(1)  J  (k,r)   -  I  Y(k,co)  exp(-ia)l)da>  =  (2n)"^  J  R(p,r )exp(ik.J)  d?  . 

The  mixed  quantity  $(k,T)  is  a  spectral  density  in  the  variable  1^,  and  an  un- 
normalized  correlation  function  in  the  variable  T;  we  shall  call  ^  (k,?:')  the 
(spatial)  semi trans form  of  the  random  process  ^  6(r,t). 

To  a  good  approximation,  the  relation  between  the  covariance  of  the 
received  scattered  current  and  the  dielectric  noise  is  given  by  the  fonrnola 

(2)  <l(t)l*(t  +  r)>  -  (nDA^RR')^  (2n)^  <;(^6)2)>  J  (",r)   , 

where  K  is  the  scattering  wave  vector,  R  is  the  distance  from  A  to  0,  and  R  is 
the  distance  from  A  to  0.  Dividing  eqiiation  (2)  by  <^ll|'^  ,  we  obtain 

(3)  c(r)  -  f  (l?,r)/f  (K,0)  . 

Eqiiation   (3)    is   the  starting  point  for  ovir  analysis  of  the  time  structure  of  l(t). 
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3,  Convection  of  the  dielectric  noise  and  its  effect  on  the  scattered  radiation 
As  just  defined,  the  dielectric  noise  ^£^  e(r,t)  and  its  semitransform 

&   (iCjT)  are  measured  in  a  coonlinate  system  fixed  vdth  respect  to  the  earth. 
On  the  other  hand,  txirbulence  theory  leads  most  directly  to  statements  about 
qi;antities  measured  in  a  Lagran£;ian  coordinate  system,  i.e.,  a  s^'stem  moving  with 
the  local  i-and  velocity  v(r,t).  We  must  now  separate  the  contribution  to 

$  (1?,  T)  due  to  convection  by  the  velocity  field  from  the  contribution  due  to 
time  variation  of  the  microstructure  of  the  dielectric  t\arbulence  as  seen  in  the 
Lagranglan  system.  Things  are  considerably  simplified  if  we  realize  from  the 
outset  that  the  convective  effect  is  due  primarily  to  the  mean  velocity  field 
and  the  macro-eddies,  and  that  the  macro-eddy  velocities  change  only  slightly 
during  the  correlation  time  of  the  received  scattered  radiation,  v/hich  experi- 
ment shows  to  be  of  the  order  of  seconds  at  most  (cf.,  e.g.,  Trolese  ^  -').  Let 
us  nov/  divide  the  scattering  volume  V  into  n  equal  cubical  subvolumes  V^,...,V 
of  side  d  satisfying  the  relation  L  «  d  <  L  ,  where  "^  L  =  2n/K.   Then  our  fading 
calculations  will  be  acc\irate  to  within  at  least  an  order  of  magnitude  if  we  re- 
place v(r,t)  in  each  siibvolume  V.  by  an  average  value  v.,  obtained  by  space- 
averaging  v(r,t)  over  V.  and  time  averaging  v(r,t)  over  an  interval  of  the  order 
of  the  correlation  time  of  the  received  scattered  radiation.  It  follows  from 
the  stationarity  properties  of  v(r,t)  that  the  random  variables  v.,  1  <  1  <  n, 
are  identically  distributed  with  probability  density  p(v)  =  p(v  ,v  ,v  ) ,  i.e., 
the  probability  that  any  of  the  v.  has  components  lying  in  the  intervals 
(v  ,v^+  dv^),  (Vy,Vy+  dVy),  iy^y\^   dv^)  is  given  by  P(\»Vy>\)-  Moreover, 
p(v)  is  approximately  the  univariate  probability  density  of  the  velocity  field 
v(r,t),  despite  the  smoothing  out  of  the  finer  fluctuations  entailed  in  the  de- 
finition of  the  v.. 

We  now  represent  C(T)  as  a  sum  of  contributions  coming  from  the  dif- 
ferent  V.,  i.e.,  we  write 

See  relevant  remarks  by  G.K.  Bate helor  W  ,  p.  120j  W.  fieisenbergL?J ,  oectlon  5j 
and  A.N.  Kolmogoroff  1^°-^ . 

The  condition  that  L  «  d  is  required  for  the  validity  of  equation  (h)  belowj 
see  [S],  In  typical  radio  scattering  experiments  the  relative  size  of  V  and 
l3  is  such  that  n  will  be  a  very  large  nimber. 

Since  the  volumes  V.  are  large  compared  to  the  size  of  the  scattering  eddies, 
there  is  negligible  correlation  betvreen  the  scattering  eddies  in  the  different 
Yj_.      Thus  in  (jj)  we  can  neglect  interactions  between  the  contributions  to  l(t) 
coming  from  different  V.. 
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(ii)         C(r)  -  (l/n)  T       I,  (K,r)/$.  (K,0)   . 
i=l   ^         ^ 
_   -* 
The  subscript  on  $.  (KjT)  indicates  that  it  pertains  to  the  subvolume  V.,  where 

the  local  xcLnd  velocity  is  approximately  v..  ^ .  (K,?')  is  to  be  calculated  in  a 
coordinate  system  fixed  mth  respect  to  the  earth.  To  separate  the  contribution 
to  \  .   (k,'?')  due  to  convection  from  that  due  to  the  time  variation  of  the  turbu- 
lence itself,  we  introduce  the  quantity  ^^{kff),   the  semitransform  of  the  di- 
electric noise  measured  in  the  moving  Lagrangian  system j  5  (k,?')  can  be  deduced 
from  turbulence  theory,  at  least  in  principle,  and  will  be  assuined  to  be  the 
same  function  in  all  the  V..   The  relation  between  ^  .  (k,T)  and  ^  (k,?')  can  be 
deduced  as  follows.  We  note  first  that  by  eq\jation  (l),  ^.  (k,T)  is  the  p-trans- 
form  of  the  correlation  function  R. (p,?^)  measured  in  the  fixed  volume  V.,  whereas 
I  (k,T)  is  the  p-transform  of  the  correlation  function  R  (p,?')  measured  in  the 
Lagrangian  system.  Let  S.  be  a  coordinate  system  moving  with  the  velocity  v. ,  and 
let  the  space  and  time  coordinates  measured  in  S.  be  r.  and  t.,  related  to  the  co- 
ordinates r  and  t  of  the  fixed  system  by 

(5)  r^  =  r  -  v^t  , 

t.  =  t       ,      l<i<n. 

The  difference  variables  p.    =     r.-  r.,     T.    =     t.-  t.   and  p=r-r,      6=t-t 

transform  in  the  same  way.     Novr  by  definition 

(6)  Ri(p,r)    =  V^i'^i^    • 

If  we  multiply  both  sides  of  (6)  by  (2n)~'^  exp(ik»p), integrate  from  -co  to  oo  in 
all  the  components  of  p  and  use  ($),  we  find 

J^(k,r)   =  (2n)"'   j  R^(p,r)exp(ik.^d?  =  (2n)'^  J  Rjp^,^  )exp(ik.^dp 

-  (2k)-^  j  Rjp.,ri)exp[ik.(p.+  v.'r)]dp.  =  exp[i(k.v.)r)  J^(k,r^). 
Thus  the  desired  relation  between  J.  (k,T)  and  $  (k,?')  is 

(7)  ii(k,r)    =   exp[i(k.v.)r]  §^(i^,r)    . 

With  the  use  of  (7),  equation  (h)   becomes 

(8)  C(r)  ~  [$^(K,r)/$^(K,0)]  (l/n)  ^   exp[i(K.v.)r]  . 
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We  assume  that  the  random  process  v(r,t),  and  consequently  the  random  variables 
v.,  have  the  ergodicity  needed  to  identify  space  and  ensemble  averages.  Thus 
the  last  factor  in  (b)  can  be  written 


(9)       (1/n)  21   e^P 
i=l 


[i(K.v.)r]  ~  I    exp[i(K-^r]p(v)  dv  =  Z(r) 


where  p(  v)  is  the  coin;non  probability  density  of  the  random  variables  v. .  The 
inte,-^ral  in  (9)  is  the  i'burier  transform  of  a  non-negative  function,  and  .lence, 
by  a  theorem  of  Boclaier-Khintchine  ,  it  is  a  correlation  function,  wliich  v;e 
have  denoted  by  %('?')•  1'hus  the  correlation  function  of  the  received  scattered 
signal  is  the  product  of  two  correlation  functions 

(ic)        c(r)  -  [j^(K,r)/ J^(K,o)]%(r)   =  ^{t,r)X(^z)  , 

the  first  of  which  represents  the  fading  due  to  time  variation  of  the  scattering 
eddies,  and  the  second  of  which  represents  the  fading  due  to  convection  of  the 
scattering  eddies  by  the  mean  wind  and  the  macro-eddies.   (The  reader  is  reminded 
that  the  product  of  two  correlation  functions  must  be  a  correlation  function  (cf. 
Loeve  i-  -,  p.  J-168).)  If  the  scattering  eddies  are  constant  in  tine,  so  that 
(^{Y.yZ)   -   1>  we  have  simply 


c(r)  ~Z(r)    -    I   exp[i(K.v)r]  p(^ 


dv 


which  is  recognized  as  the  standard  expression  for  the  Doppler  fading  of  signals 
retvirned  by  raiidomly  mo'idng  aggregates  of  discrete  scatterers*- -*' ^-^ . 

To  derive  expressions  suitable  for  comparison  with  experiment,  we  pass 
from  the  complex  random  process  I(t)  to  its  absolute  -/alue  or  envelope  |l(t)I. 
Since  l(t)  is  a  complex  gaussian  process,  with  identically  distributed  and 
orthogonal  real  and  imaginary  parts'  ,  it  follows  that  |l(t)|  is  a  Rayleigh  process. 
3y  this  we  mean  that  the  distribution  of  the  random  variables 

|i(t^)|,  ...  ,  iKtJI 

for  any  t^f.,.ft   .   and  n,  is  the  appropriate  miiltivariate  Rayleigh  distribution 
given  by  Hoffman  ■=■  -^ .   The  envelope  correlation  function  is  the  quantity 

c^ir)   -   <(!i(t)|  -  <|il>)(li(t^r)i  -  <|ii»>    . 


■jrtf-, 


^See  M.  Lo^veM,  pp.  207,  U82  and  B.V.  Gnedenko '-"'■^J  . 

These  properties  of  the  real  and  imaginary  parts  of  l(t)  are  a  consequence 
of  the  reality  and  staticnarity  of  ^e(r,t)  . 


It  is  easily  shown'-  -I  that  the  relation  between  C  (T)  and  CCT)  is 


(II)       C  (T) 


[Jc'('?')-ic^('^)*-] 


We  now  choose  as  a  convenient  measiire  of  the  envelope  fading  rate  the 
quantity  Z,  the  expected  number  of  zero  crossin^per  second  of  the  centered  random 
process    Il(t)  |   -    <(|ll>     •     %  a  well-knoim  formula  L-   J 

Z     =     (1/n)     /-C^'  (0)'    , 

where  the  prime  denotes  differentiation  with  respect  to  ?'.  By  (lO)  and  (11) 

(12)  Z  ^r    -/-/   (K,0)  -  %"(0)'  ,     Y  =   [(2n)a,-n)]-l/^   , 

where  we  have  used 

0   (K,0)  '     X{0)     =1    (normalization), 

/(K,0)  =  ;;^:'(0)   =  0    (oddness)  . 

To  evaluate  %    (o),  we  expand  the  exponential  in  the  defining  integral  (9)  and 
find  that 

(13)  -  X  (0)   -   f  (K.^S(^  dv  -^<(K.t)^>    . 

In  radio  scattering  experiments  K  is  usually  directed  vertically  doimward  (see 
Figure  1).  Thus,  combining  (12)  and  (13),  we  have  finally 

(111)  Z  -  Y  /-12f"(K,0)  +<v^>  K^'   . 

If  the  scattering  eddies  do  not  vary  in  time,  equation  (lU)  reduces  to 

(1^       Z  --  Y  <\y^^   K  ~  (2n©Y/c)  <v^>^/^  f  , 

the  familiar  linear  Z-versus-f  law  characteristic  of  Doppler  sliifting.  Moreover, 

,11  ..*  . 
if  the  time  variation  of  the  scattering  eddies  is  such  that  $?  (K,0)  is  propor- 

2 
tional  to  K  ,  i.e., 

-  S2f"  (K,0)  =  a  K^  , 


then   (lii)   becomes 

(16)  Z    -   (2n&Y/c)  yfT<7>~'  f    . 

It  is  important  to  note  that  (l5)  and  (16)  are  not  only  the  eq\iations  of  straight 
lines,  but  of  straight  lines  passing  through  the  origin  of  the  (Z,f)  coordinate 
system.  Equations  (l5)  and  (l6)  correspond  to  the  resiilts  of  earlier  work  l- -' ' '- -^ , 
where  the  dielectric  turbulence  itself  was  represented  heuristically  as  a  cloud 
of  independent  discrete  scatterert.  As  we  shall  see  below,  the  tiorbulent  mixing 
model  of  dielectric  turbulence  leads  to  a  different  K-dependence  for  ^  (K,C)  and 
to  a  different  Z- versus -f  law, 

ij.  Fading  rate  with  turbiilent  mixing  model  of  dielectric  turbulence 

To  apply  statistical  turbulence  theory  to  an  analysis  of  ^  (K,0),  we 
shall  assume  that  the  dielectric  noise  is  produced  by  turbulent  mixing  of  dielectric 
constant  inhomogeneities  .  Guided  by  the  work  of  Obukhof  f '-  -^  and  Corrsin  l-  -' ,  we 
assume  that  the  turbulently  mixed  field  ^  e(r,t)  has  the  same  mic restructure  as 
the  underlying  turbulent  velocity  field  v(r,t).  VJith  such  a  model  we  are  justified 
in  inferring  the  time  structure  of  dielectric  constant  eddies  of  a  given  size  from 
the  time  structxu'e  of  velocity  eddies  of  the  same  size.  On  the  other  hand,  if 
dielectric  turbulence  were  due  to  some  other  mechanism,  e.g.,  the  density  fluctua- 
tions associated  with  turbulent  pressure  fluctuations,  we  could  not  make  this 
identification.  Indeed,  in  the  case  of  the  pressure  fluctioaticn  model,  now  abandoned 
as  an  important  source  o£   dielectric  turbulence,  the  time  structure  of  pressure 
eddies  of  a  given  size  is  related  in  a  complicatea  way  to  the  time  struct\ire  of 
velocity  eddies  of  many  different  sizes  ^  -'. 

Let  J  (kjTlO  be  the  semi  transform  of  the  turbulent  velocity  field  de- 
fined by  (1),  with  T(k,o>)  and  R(p,T)  interpreted  as  the  spectral  density  and  corre- 
lation function  of  the  velocity  field.  Despite  the  fact  that  a  satisfactorj'  deriva- 
tion of  J  (k,T)  has  not  yet  been  achieved  (except  for  T  =  O),  we  can  infer  the 
quantity  J  (k,0),  and  presumably  0     (k,0),  to  within  a  constant  of  proportionality, 
by  using  dimensionality  and  similarity  considerations.  The  argument  is  familiar 
from  the  work  of  von  Weizsacker  I-  -^  and  Heisenberg^'^^-',  and  goes  as  follows.  If  the 

Reynolds  number  of  the  turbulent  flow  is  large,  the  rate  at  which  kinetic  energy 
■{t-  — ^— — — ^— ^^— ^— __^— ^_^__^_— _^— — -^— _— — _ 

The  raicrostructure  of  dielectric  turbulence  is  still  a  veiy  open  issue.  Thus, 
■the  considerations  of  this  section  are  perforce  of  a  speculative  nature.  Evi- 
dence supportinij  the  view  that  dielectric  turbulence  is  due  primarily  to  turbu- 
lent mixing  of  djTiamically  passive ^quantities  like  temperature,  humidity,  and 
electron  density,  is  presented  in 


^ 


-  10  - 


is  transferred  f^om  eddies  larger  than  A  "   2n/k  to  eddies  smaller  than  -C  is 
independent  of  J?  ,  provided  that  J?  lies  in  the  inertial  range  L' -1  >  L  J ,  Thus,  if 
Vg     is  the  velocity  of  eddies  of  size  ^  ,  we  have 

(17)  "^  /^   "  ''°'^^' 

Defining  a  characteristic  time  2^    =    A/v      for  the  eddies  of  size  £  ,  and  using 
(17),  we  find 

-^  -.  const,   i^/^  . 

We  now  assume,  following  Heisenberg,  that  ^  (k,  j?  )  is  a  universal  furiction  of 
the  dimensionless  variable  T/Zj    ,   which  is  proportional  to  k  '    t»     (The  fact 
that  J  (k,T)  does  not  depend  on  the  direction  of  k  is  a  consequence  of  the 
local  isotropy  of  eddies  in  the  inertial  range  l- -^ '  ^ -'. )  Thus,  assuming  that  the 
time  structure  of  the  dielectric  constant  eddies  'parallels'  that  of  the  velocity 
eddies,  we  can  write 

The  desired  expression  for  ^  (K,0)  is  therefore 

(18)  f   (K,0)   =  g"  (0)  K^/^  . 

It 

Since  correlation  functions  are  concave  downwcird,  we  know  that  g  (o)  is  negative, 

but  the  exact  value  of  g  (O)  cannot  be  determined  without  a  more  quantitative  theory 
If  we  substitute  (18)  in  (Ih),  we  obtain 

(19)  z  ~  r  /-g"  (0)  k"/''  *  Ol>'  ' . 
J'    [-.■'(0,/<.^>j3A, 

we  can  write  (19)  as 

(20)  Z  ~  Y  Ol>^^^     ^o    /(KAo^^'^^^  (kAq)^'  . 

The  validity  of  (20)  depends  on  two  conditions,  which  are  equivalent  to  the  re- 
quirement that  the  scattering  eddies  lie  in  the  inertial  tange.  The  first  is 


Defining  a  wave  number  K  by 
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that  L  =  2n/i(  ~  \/Q   satisfy  the  inequality  L  «  L  ,  i.e.,  that  the  size  of  the 
scattering  eddies  be  mijch  less  than  the  macro-eddy  dimension.  This  condition 
seems  to  be  met  in  raost  radio  scattering  experiments,  both  in  the  ionosphere 
and  in  the  troposphere.  The  second  condition  is  that  L  »  "^  ,  where  V   is  the 
characteristic  dimension  of  the  smallest  eddies,  which  is  determined  by  the 
action  of  molecular  dissipation  processes.  This  condition  does  not  seem  to  be 
met  in  most  ionospheric  scattering  experiments".  Thus,  although  the  results  of 
previous  sections  apply  to  both  ionospheric  and  tropospheric  scattering,  the  ap- 
plication of  equations  (17)  through  (20)  is  limited  to  tropospheric  experiments, 
where  the  condition  L  »  ''^  seems  always  to  be  satisfied. 

An  alternate  form  of  (20)  in  terms  of  the  carrier  frequency  is 

(12)       Z  -(2n©Y  <0^/Vc)  f;)/(f/f^)^/^  ^  (f/f^)''  , 

where  now 

f  =  (c/2ng)  K  , 
0  o 

In  figure  2  we  plot  the  function  h(x)    =     )oc ''    +  x  ,  wMch  is  the  Z-versus-f  law 

if  Z  is  neasiured  in  units   of   (2n©Y  <s,^    /*      /c)f     and  f  is  measvired  in  units  of 

f   ,     Noticeable  departures  from  the  Doppler  laws   given  by  eqtiations  (l5)   and   (l6) 

o 
can  occur  in  two  ways.  First  of  all,  if  (21)  is  valid  at  frequencies  f  ^  f  /2, 

which  it  may  not  be   ,  then  at  such  frequencies  the  Z-versus-f  law  is  nonlinear 
(see  part  of  graph  near  origin).  On  the  other  hand,  even  if  (21)  does  not  apply 
at  frequencies  f  <,  f  /2,  so  that  (21)  is  valid  only  for  (Z,f)  points  lying  on 
the  portion  of  h(x)  wliich  is  closely  approximated  by  a  strait^ht  line,  the  con- 
tinuation of  this  straight  line  vri.ll  not  pass  through  the  origin,  as  predicted  by 

See,  for  example,  F.  Villars  and  V,F.  Weisskopf  L^PJ  and  A.D,  WheelonL^lJ. 

For  estimates  of  V   in  the  troposphere,  see  A,M,  Cbukhoff '-  -•  and  G,K. 
Batchelor  ^^-^ , 

Corresponding  to  the  condition  L  «  L  or  KL  »  2n,  equation  (21)  is  valid 
only  for  frequencies  such  that  OfL  /c  »  1.  Thus,  if  a  noticeable  departure 
from  linearity  is  to  be  predicted  by  (21),  we  must  have  ef  L  /c  »  1.  Whether 
or  not  this  condition  is  met,  a  crude  estimate  shows  that  f^  is  a  frequency 
of  the  order  of  megacycles,  (i.e.,  comparable  to  radio  scattering  frequencies) 
and  that  the  fading  rate  Z  is  a  frequency  of  the  order  of  several  cycles  per 
second  (i.e.,  comparable  to  observed  fading  rates). 
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the  Doppler  laws  (l^  and  (l6)  (see  broken  line  in  Figure  2).  Thus,  in  any  case, 
time  variation  of  the  scattering  eddies  shoiild  produce  a  measurable  effect  on 
the  time  structure  of  radio  waves  scattered  by  dielectric  turbulence  in  the  tro- 
posphere • 


h(x) 


2.0 


0.5- 


Figur 


e  2 


Graph  of  the  function  h(x)  «    J  iT'^  *    x     . 


13  - 


5,  Conclusions 

We  have  shown  that  the  time  structure  of  electromagnetic  radiation 
scattered  by  a  continuous  dielectric  constant  field  contains  contributions  from 
two  sources,  vne  is  the  time  variation  of  the  spectral  density  of  the  random 
process  at  the  scattering  wave  number,  as  seen  in  a  system  moving  with  the  local 
wind  velocity,  which  in  the  case  of  atmospheric  turbvilence  is  itself  a  random 
process.  The  other  is  Doppler  shifting  due  to  convection  of  the  dielectric  con- 
stajit  patterns  that  cause  scattering,  the  so-called  scattering  eddies,  by  the 
mean  wind  and  by  the  large  velocity  eddies.  Our  result  is  summarized  in  equation 
(10),  which  represents  the  correlation  function  of  the  received  scattered  cur- 
rent as  the  product  of  two  correlation  functions,  one  representing  time  variation 
of  the  scattered  eddies,  and  the  other  convective  Doppler  shifting. 

On  less  certain  grovmd,  we  have  assumed  that  dielectric  turbulence  is 
produced  by  the  turbulent  mixing  of  dielectric  constant  inhoraogeneities,  and  obeys 
the  Kolmogoroff-Obukhoff  structure  law.  Although  a  complete  solution  of  the  space- 
time  structure  of  turbulence,  in  particxilar  of  dielectric  turbulence,  is  not  yet 
available,  it  is  possible,  by  the  use  of  familiar  dimensionality  ajid  similarity 
arguments,  to  deduce  the  dependence  of  the  fading  rate  of  the  received  envelope 
on  the  carrier  frequency.  The  result  is  sunmaiuzed  in  equation  (21)  and  figure  2, 
which  apply  only  to  tropospheric  scattering  and  which  exhibit  significant  depar- 
t\ires  from  the  Doppler  law  that  goes  with  the  scattering  blob  models.  Experimental 
evidence  in  favor  of  the  fading  rate  law  of  (2l)j  as  opposed  to  the  Doppler  law, 
appears  indecisive.  For  example,  Trolese  •-  -'  reports  simultaneous  measurements 
of  the  fading  rate  (defined  somewhat  differently)  at  three  frequencies,  the  minimum 
number  of  points  needed  to  decide  for  or  against  a  linear  law.  He  notes  that  of 
ills  two  sets  of  data,  consisting  of  three  (Z,f)  points  each,  only  one  set  fits  fairly 
well  on  a  straight  line  through  the  origin.  It  is  clear  that  a  decisive  test  of 
(21)  awaits  more  extensive  experiments,  which  it  is  hoped  this  paper  may  in  some 
measure  stimulate. 
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Abstract 

The  fading  of  radio  waves  scattered  by  dielectric  tiirbvilence  is 
shown  to  be  the  result  of  two  effects.  The  first  is  time  variation  of  the 
scattering  eddies  as  seen  in  a  coordinate  system  moving  with  the  local  vrind 
velocity.  The  second  is  Doppler  shifting  prodxiced  by  the  convection  of  the 
scattering  eddies  by  the  mean  wind  and  by  the  macro-eddies.  In  the  tropo- 
sphere, the  scattering  eddies  lie  in  the  inertial  range  of  statistical  tur- 
bulence theory.  This  makes  it  possible  to  find  the  envelope  fading  rate  of 
the  received  scattered  signal  to  within  a  constant  of  proportionality  by 
using  dimensionality  and  similarity  arguments.  The  result  is  an  expression 
for  the  fading  rate  as  a  function  of  carrier  frequency  which  departs  signifi- 
cantly from  that  found  with  'scattering  blob'  models. 
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